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ON THE IRREGULAR INTEGRALS OF LINEAR DIFFERENTIAL 

EQUATIONS.* 

By Clyde E. Love. 

1. In the present paper it is proposed to study the integrals of a homo- 
geneous linear differential equation in the neighborhood of an irregular 
point, particularly with a view to deterrnining the existence and form of 
their asymptotic developments.! The existence of such developments is 
well-established provided the roots of the so-called characteristic equation 
are all distinct,! but for the case of multiple roots no general discussion has 
appeared. In the following work, while the roots are unrestricted as to 
their order of multiplicity, a certain limitation is placed upon them for the 
purpose of reducing the algebraic difficulties of the analysis. This limita- 
tion is of such a nature that the case of distinct roots is included as a special 
case of the problem here treated. 

The irregular point is taken as the point at infinity, and the independent 
variable is restricted to real and positive values. 

2. As a prehminary step we shall formulate two general theorems§ on 
linear differential equations, which will form the basis for the subsequent 
investigation.. 

Take for consideration the differential equation 

(1) y M + a 1 (x)y^~» + ■ ■ ■ + a n (x)y = 0, 

whose coefficients a,i(x), a 2 (x), •••, a„(x) are continuous and possess n 
continuous derivatives for all sufficiently large positive values of x. Let 
Zi, z 2 , • • •, z„ denote n auxiliary functions of x which, together with their 
first n derivatives, are likewise continuous when x is large and positive, 
and are such that the determinant 







zi zi ••• Zi (B-1) 


(2) 


Q(*) = 


Z 2 Zi ••• z 2 (n-1) 

y v ' ... v (n— 1) 
An An An 


never vanishes. 


Let A r (x) 


be the minor of Q(x) 



* Read before the American Mathematical Society, Dec. 27, 1913. 
f In PoincarS's sense: cf. Acta Mathematica, vol. 8 (1886), p. 297. 
X Cf., for example, Horn, Journal fur Mathematik, vol. 138 (1910), pp. 159-191. 
§ Cf. Dini, Annali di Matematica, Ser. Ill, vol. 2 (1899), pp. 297-324. 
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2 r (B-1) . Also define 

Z r (x) = Zrdn - (*#+.{)' +•••+(- 1) "^{z^) <-» + (-l)"S r <»>, 

(3) r- 1,2, ••.,!», 

and denote by q(x, x x ) the determinant formed from Q(x) by replacing 
the element 3 r (n_1) (x) by Z r (xi), r = 1, ■•-, n. 
Place 

(4) X(x, xi) q^ 

(— l) n-1 C 4 r (i) 

(5) 9r(x) = U T , o(x) = V r , o(x) = Q/^ , r = 1, • • • 71, 

where C r is an arbitrary constant; 

••• I K(x, x 1 )K(xi, x 2 ) 

• • • X(x A _„ x x )sr r (x A )dx x • • • dx 2 dxi, X = 1, 2, • • • , 

Vr, *(X) = ff-f £(*, Xi)K(Xi, Xt) 

• • • X(x x _„ x A )gr r (x x )dx x • • • dxidxi, X = 1, 2, • • •. 

Then we have 

Theorem A: Suppose that for all values of x greater than some constant, 
the series 

00 

(6) 2/r(x) = Z V r , a(x) 

satisfies the following conditions: (a) the series converges; (6) the series for 
y r (xi) when multiplied by K{x, Xi) may be integrated term by term with respect 
to Xi from x to °o ; (c) i/ie series defines a function y r (x) such that each of the 
integrals 

/»00 

J y r (x)Z ,(x)dx, s = 1, • ■ ■, n, 

/ias a meaning. Then for such values of x the function y T (x) is an integral of 

(I)- 

Theorem B : Suppose that a constant a can be found such that for all 

values of x > a the series 

00 

Vr(x) = Z«r, x(x) 

satisfies the following conditions: (a) tfie series converges; (b) the series for 
y T (xi) when multiplied by K(x, Xi) may be integrated term by term with respect 
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to X\ from a to x. Then for such values of x the function y r (x) is an integral 
of (1). 

To prove Theorem A, place 



(7) 



p., o = «., 

P», 1 = 2«ai — 2»' = 2«Ol — P». </, 

■ P», n-1 = Z.fln-1 — P., n-t'', 3 = 1, 



■,n; 



2/rZ.dx, s = l, • • •, r - 1, r + 1, 
y r Z r dx + C r ; 

A r (x) = 



n; 



Pi. Pi, 1 • • 

Pn. Pn. 1 ' * 
Pi. Pi. 1 • 
Pn. Pn, 1 * 

Now by condition (&) equation (6) may be written 



A(x) = 



Pl, n-2 


#i(x) 


Pn. n-2 *»(Z) 


Pl. n-1 


= (" 


Pn, n-1 





»(n-l) 



(8) 



Xeo 
y r (xi)K(x, Xijdxt. 



By substituting the values of gr r (x) and K(x, Xi) in (8) we find 

y r (x) = A r (x)/A(x), 

so that it suffices for our proof to show that this function is an integral of (1). 
To do this, consider the system of n functions 770, 171, •••, »?n-i each 
defined for all values of x sufficiently large by means of the following 
system of n linear equations: 

(9) p., oln-l + P». 1 Vn-2 + • • • + P,, n-1 Vo = **(z), S - 1, • • • , TO. 

We have at once 

A r (x) 



(10) 



570 = 



A(x) 



= 2/r 



Upon differentiating equations (9) with respect to x, and making use of 
(7) and (10), we find 



(11) 



Z.0 + P>, </01 + P», l'02 + • • • + P., n-i'dn-l = 0, 8=1, 



•, TO, 
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where 

(12) B = j?„_i' + ai(x)vn-2 + •■' + a n _i(z) V + a„(x)r) , 

Si = IJn— 1 — Vn—2 t 

(13) • #2 = Vn-2 — Vn-s', 

■ 0n-l = Vl ~ Vo • 

The system (11) consists of n homogeneous linear equations in the n 
quantities 6, 0i, • • •, 0„_i. By virtue of the relation 

Pt, ,' = z«a«H-i — Pt. s+i, t = 1, • • •> n; s = 0, 1, • • •, n — 2, 

the discriminant of the system reduces at once to (— l) B-I Q(a:), and hence 
by our hypotheses does not vanish for any value of x under consideration. 

Whence 

8 = 0i = ■ • • = 0„_i = 0, 
or by (10) and (13), 

V, = j/r 00 , s = 0, 1, ---.n. 

Substituting in (12), we find 

y r ^ + atWyS"-" + ••■ + a n (x)y r = 0, 

which was to be proved. 

The proof of Theorem B follows similar lines, and may be omitted. 

These theorems differ only slightly from certain results due to Dini.* 
He was first led by synthetic processes to the converse theorems, and was 
then able to establish the direct theorems by reversing his previous reason- 
ing. The appearance of artificiality in the proof as here outlined is chiefly 
due to the fact that only this reversed line of argument is presented. 

3. In the differential equation (1), suppose that the coefficients ai(x), 
02(x), • • •, a„(x) are real or complex functions developable, for large real 
positive values of x, in asymptotic (or convergent) power series of the form 

c r (x) « x Tk c r , o + -^ + • • • i t = 1, • • •, n, 

k being or a positive integer. Suppose also that the first n — r derivatives 
of a T (x) possess asymptotic developments in the same z-region. 

We proceed to apply the theorems of § 2 to the problem of fi n di n g 
asymptotic solutions of (1), valid for large real values of sr. 

Suppose that the characteristic equation 

<f>(m) = m n + d, ow" -1 + • • • + a„, o = 

corresponding to (1) has I different roots m u m 2 , • • • , mi, occurring n u m, 
' • ', ni times respectively: (»i + n 2 + • • • + rij = n). 

* Loc. cit. 
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Let us choose n auxiliary functions (cf . § 2) of the form 

(14) z» 1+ ... + »,. 1+s - x-"-. «e-*,*»**..<w r = l, . . . , I; q _ 1, . . . , n,, 

where 

1 2 J. 

/'••(*)- TTFT+—; — r + rT7~^ + '" +_ i~ 



and 



*r.,(s)-- t f- + — f— + ••• + -^rr' 



<*>-.«> Pr.q,- • •) v r,q being undetermined constants and s an arbitrary positive 
integer. 

For brevity let us write 

2 »i+ •••+»,_,+« = Zr,q, 

Zni+'"+ n r-i+Q ~ "r,qi 

?«,+•• •+n r ^ 1 +q — Qt, «• 

Now upon forming Z r , 9 (x) as given by (3), we find that 
(15) Z r ,q(x) - (- l)"* r .,a;"*[S r „(*) + ar-Vi,,.,^/*)], 

where* <S r , g (x) is a polynomial of degree n r (s + A;) — 1 in x~ 1 ' nr . In 

<S r , q(x), the coefficients of ar''/*-, j = 0, 1, • • •, n, — 1, vanish identically 
because of the fact that 

<p(m T ) = <p'(m T ) = ••• = ?><»•-" (m,) = 0. 

Further, writing this polynomial in the form 

S,.,(x) = ar^A,,, + Br,,*- 1 /"- + r P . «*-»/*• + • • • + N,,,ar*+«/~], 

we have 

g<->Qn,) > . 

Ar,« = — "i <*,,,*' + f{m r ), 

where 

yp{m) = Ci, im n_1 + a 2 , iTO" -2 + • • • + a«, ij 

* Throughout the work, the symbol P(x) will denote a function expressible, for large values 
of x, in the form 

P(s).«.+g + ... + a ' + J' (X \ lim «,(*)-<), 
where p is an arbitrary positive integer. 
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and 



<p (nr) (ro r ) 
(n r -l)!' 



B r,« = 7^— ZTfi ar . «"' ^r.a + "r,,(fflr, «r, «), 



IV, 9 = ( w _ 1 \j «r,g" r ~ 1 7r, g + T r , q (m r , OL r ,q, fir.q), 

Nr '« = 7^ TTi"''-'"' -1 ''' - -'' + "r.aC 771 "-. a r,q, * * ", Mr, j), 

where <r r , ,, r r , «, • ■ •, « r , s are certain polynomials in the indicated argu- 
ments, which it would be useless to write out in full. 

Let us now try to determine the coefficients a r , q , fir.q, •••, v T , q by 
placing* 
(16) " A,, , = B r , ,= •••= N,. , = 0. 

For simplicity we shall assume that no multiple root of the characteristic 
equation is also a root of the equation \p(m) = 0, or in other words that there 
is no value of m satisfying simultaneously the three equations! 

(p(m) = <p'(m) = \p{m) = 0. 

Under this hypothesis equations (16) can always be solved, and they 
serve to determine n definite, distinct functions z r , 9 of the form (14). 

4. The quantities R[f r , q (x)], where R[x] denotes the real part of x, 
will play an important part in what follows. We consider first the case 
in which these real parts are all distinct, and shall suppose for denniteness 
that when x is sufficiently large 

R[fi, i(x)] < R[fi, .(*)] < < R[fi. »,(*)] 
(1?) < R[h. i(x)] < < R[fi. .,(*)]. 

The argument to be used will apply with changes only in notation if 
these real parts occur in any other order when arranged with respect to 
their magnitude, so that the assumption (17) involves no loss of generality. 

Let n' denote the least common multiple of n u n 2 , • • • , n h and set 

fix) - /l, l(x) + fl, .(*) + • • • + fl, n,(x), 

(n — l)nk — n + I 
r\ = 171, i + iji, 2 + * * " ~t~ vi, n, — 2 ' 

A = (- 1) 2 n(m P - w«) n,B, n(a r , 9 - a r , p ), 



* As soon as the constants v r „ have been determined, s is to be selected so as to satisfy the 
inequality (21). 

t It is understood, of course, that the case <p(m) m is also excluded, except when the roots of 
the characteristic equation are all distinct. 
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where the first product is composed of all the factors m r — m 3 that can be 
formed with r = 2, 3, ■ • • I and r > s; the second consists of all factors 
a r , g — a,, p that can be formed with r = 1, 2, • • •, I; q = 2, 3, • • ■, n T 
and q > p. 

With this notation we find, by (2), 

Q(x) = Ae-^x-'PaCx 1 "*'), 

where the constant term in* P 2 (x 1/n ') is unity, so that the hypothesis 
Q(x) =t= for all values of x under consideration is satisfied. Further, with 
proper choice of the arbitrary constant C r , g , we may write by (5), 

(18) g r . q (x) = e'r, ,<»>afr. «P,. ,, ,(**'"'), 
where we have set 

Pr, q = Vr, q + 1 — (n — l)k, 

and where the constant term in P 3 , r , 9 (x 1/n ') is unity. And by (4) we have 

(19) K(x, xi) = fi fie''. 9 (l) x"'. *P 3 , r, ,(x 1 /»')e-''- « (x,) xr"'. «-'P 4 , r. ^Xi 1 /"')- 

5. We shall now write out by Theorem A a trial solution of (1) corre- 
sponding to that one of the functions g r , 9 (x) whose exponential factor has 
the least real part, which by (17) is g u i(x), and show that the conditions 
of the theorem are satisfied. We have 

(20) t/i, i = gi, i(x) + £ Vi, i, x (x), 
where 

»1. 1, x(x) =1 I • • • I •£(«, Xi)Z(Xi, x 2 ) 

• • • X(xw, xOffi. i(x0cte x • • • dx 2 dxi. 
Placing 

?r, «(x) = e'r. q^X'T, 9, 

I n r 
K(X, Xi) = 2 Tle'r. q^Xfr, qe~ f r, qWxr'r. «~», 

Jf»o0 /*OD /*OD __ 

I • • • I K{x, Xi)K(x lf x 2 ) 

• • • K(x K - u x K )g u i(x K )dx K ■ ■ • dxidx u 

let us study the behavior, for large values of x, of the function Vi, i, x (x) . 
We first choose an arbitrary positive integer p, and then take s so large 

* Cf. footnote, p. 149. 
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that* 

(21) s > v + 2 + 2p, 

where p is a positive number at least as great as any of the quantities 

\R[Pr, ,]|. 

Since, by (17), 

R[fi.M ~fr, ,(*0] ^0 

for all values of r and q, each of the functions 

\e~ f r, t<- x Jx~''r, g~'e f i< i^x/i. i| 

decreases monotonically for all values of x A under consideration, and we have 

\K(x^.i, xjgi, i(x K ) | < n | eA. if^-i)^-!"!. i|xr p-2 . 
Similarly, 

|X(x A _ s , x A _!)Z(x A _i, x A )ffi, i(x A )| < n 2 |eA ) i(^-2»x A _ 2 ' , ia!x A _i-' , - 2 x A -p- 2 . 
Proceeding in this way, we find 

| ••• I (x!X 2 •••x A )-"- 2 dx A 

(22) Jxi J '*~ l nA j 

• • • dx 2 dxi < |eA. i<*Vi. i| (p * 1)A • x r^^ • 

By virtue of the defining property of the functions P(x) (see footnote, 
p. 149), we may take x so large that 

|«i. i.»(*)l < 2|«i, i, A (x)|. 

00 

Whence, by (22) the series £>i, i, A (x) converges absolutely and uniformly 

for all sufficiently large values of x, so that condition (a) of Theorem A is 
satisfied. Further, by comparison with the geometric series whose general 
term is twice the right member of (22), it appears that 

2n 



Z I pi. i. *(*) i < I e'i.i<*>x>i.i| — L .-^ 



n 



(23) (P + 1 ) x * + ' 1 

■ , Ml i 4n 1 
< | c / 1>l(lVl!l | ___._. 

That condition (b) is satisfied follows from the results just established, 
together with the following theorem :f 

* Cf. footnote, p. 150. 

t Bromwieh, Infinite Series, p. 453. 
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If 2/„(x) converges uniformly in any fixed interval a ^ x ^ 6, where b 
is arbitrary, and if (p(x) is continuous for all finite values of x, then 

<p(x)2f„(x)dx = 2 I <p(x)f n (x)dx, 

|p(x)|2|/„(x)|dx is convergent. 
For, if we take a = x, x = x h <p(x{) = K(x, x{), /»(xi) = v lt i, A (xi), 

/•OO OB 

it only remains to show that I \K(x, x0 1 £ | »i. 1. *(xi) |dxi converges. But 

Jx \=0 

by (4) and (22) this breaks up into n integrals whose convergence is assured 
by that of the integrals I |e'i,i( Xl > _ As 9 (I1) :r 1 ' > i.i - ' > r. <r'|dxi, and each of 

Jx 

these converges by reason of (17) and (21). 

These same inequalities are sufficient to show that condition (c) holds, 
if it be recalled that, in (15), we have made 

S T , ,(X) = 0. 

Hence y\, i as given by (20) is an actual solution of (1). 

6. When we attempt to write out a particular integral of (1) corre- 
sponding to any other of the functions g r , 4 (x), the above process fails for 
the reason that R[f r , «(x) — f lt i(x)] > and some of the integrals in 
v r, q. k{x) diverge. To avoid this difficulty we shall first modify slightly 
our system of auxiliary functions and then use Theorem B. 

The adjoint equation corresponding to (1), viz. 

(24) z (n > - (a 1 z)<"- 1) +••• + (- l) n za„ = 0, 

satisfies all the conditions that have been imposed upon (1) ; it has therefore 
a particular integral given by (20). The functions /,, «,(x) for (24) are 
obviously the same as for (1), except that all signs are changed. The one 
having the least real part is therefore — fu n,(x), and the corresponding 
solution of (24) is found by direct computation to be 

zi, n«x = e~'i. »<<*>a;-". "<P S , j, »(xw»). 

In our system of auxiliary functions let us replace s t , », by this function 
zi, »„ thus making Z t , „, = 0. This leaves Q(x), g r , q (x) and g r , ,(x) un- 
changed in essentials, while in it~(x, Xi) and K(x, Xi) the term for which 
r = I and q = rii no longer appears. Now write by Theorem B a trial 
integral of (1) of the form 

CO 

(25) yi, n , = gi, n/x) + Z Ui, „„ A (x). 
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Put 

I • • • I K(x, xJKixx, x 2 ) ■ ■ K (x A _i, xjg r , c (x x )dx A - • • dx^dx^. 
Now we have 

R[fl. n t (X.) - fr, t (x,)] > 

for all values of r and g occurring in K(x^ u x K ). Thus each of the functions 
\ e / i,n{ x Ki-Sr,<{*>.) Z} t>i,nr<>r,<,-*] i increases monotonically throughout the interval 
a ^ x K ^ x A _i, so that 

\K(x k - h xjgt, n ,(a; x ) | < n le^^-^./'^r^+v \ . 
Similarly, 

\K(x K -2, x^Kix^, x x )gi, ni (x K )\ < n 2 |e'''"' ( **- 2 Vi p <-»- 2(p+2> I • 
And 

I u r , q , K (x) | < n A leAn^x"'.^ \x-^+» f [••■ f K1 dxi ■ • • dxidxx 

Whence by reasoning similar to that used in the previous case it appears 
that the conditions of the theorem are satisfied, and yi, „, is a solution of (1). 
7. We may now form by (25) a second solution of (24), which takes the 
form 

«l.l = C- / i,i«I-'i,iPs, 1 , 1 (l 1 '" 1 ). 

This being introduced^ into the auxiliary system in place of Z\, i, the 
term in K(x, Xi) and K(x, Xi) for which r = 1 and q = 1 disappears, since 
now Zi, i = 0. 

Upon writing out by Theorem A the function* 

00 

Vi.t = gi.i(x) + X>i,2,a(x), 

the difficulty which formerly presented itself is seen to be obviated by our 
more favorable choice of auxiliary functions, since for all values of r and q 
now occurring 

R[fl. 2(X) - fr, q (x)] S 0. 

Having proved by the reasoning of § 5 that y lt 2 is an integral of (1), 
we need only introduce a solution zi, nr -i of (24) in place of zi, n ,-i to be 
able to obtain by Theorem B another solution of (1). 

Continuing in this way, we get n integrals of (1) having one or the other 

* Of course if n t = 1, we should here form y t ,i. 
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of the forms 



(26) 



= Qt 



Vr. « = gr, q (x) + T,V T , g , A (x). 



It is easily shown, and is in fact almost self-evident, that these integrals 
are linearly independent. 

8. In case the quantities R[f r , «(x)], while all distinct, occur in another 
order than that of (17), the above argument evidently serves without modi- 
fication, merely giving us the integrals in a different order. 

In case the real parts are not all distinct, we can obtain at once by 
Theorem A particular integrals corresponding to each of the functions 
St, q (x) whose real parts are equal and least. A group of solutions of (24) 
can then be formed, after which Theorem B serves to determine a group of 
solutions of (1) corresponding to those functions f r , q (x) whose real parts 
are equal and greatest. The rest of the process is obvious. 

In particular, if the quantities R[f r , «(#)] are all equal, the complete 
system of integrals can be written down by Theorem A, at once. 

9. It is an immediate consequence of (23) and the similar result for 
the other solutions (26) that 

]hnxP \ yr ' q a~M q(X) ] = > r=l,...,l;q = l,..;n r . 
Whence 

Vr, <, ™ g r . ,(«). 

By direct substitution in the differential equation it appears that we may 
write 

where the A r , q , ,, ,-, are determinate constants and A r , «, o, o = 1- 

10. In summary, we may state the following result: 
Suppose that the differential equation 

(1) y in) + a t (x)y< n -» + • • • + a n (x)y = 

satisfies the following conditions for all sufficiently large real positive values ofx: 
(a) The coefficients ai(x), •••, a n (x) and their first n derivatives are 

continuous. 

(6) The same coefficients are developable in asymptotic (or convergent) series 

of the form 

a r (x) « x rk I a r , o + ~" H J , r = 1, • • •, n, 
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where k is or a positive integer, and the first n — r derivatives of a r (x) also 
possess asymptotic developments. 
(c) The characteristic equation 

m n + a u m n ~ l + • • • + a„, = 

has I roots mi, • • • , m i, occurring respectively n t , • • • , ni times : (tti -\ \-%i = n) 

and such that no multiple root of the characteristic equation is also a root of the 
equation 

fli. \m n ~ x + a 2 , iTO" -2 + • • • + a„, i = 0. 

Then equation (1) possesses a fundamental system of solutions y r , , (r = 1, 
' ' ' > l> Q — ]■> ' "t n r) developable asymptotically in the form 



e f r, «<*> X) X"r, «->' Z — ; 



■r, g i, j 



where f r , , (x) is a certain polynomial of degree n r (k + 1) in x llnr , the quantities 
p, t , and A T , q , », y are determinate constants, and A r , 8 , o, o = 1. 
TJnivebsitt of Michigan. 



